Insulating behavior with spin and charge order in the ionic Hubbard model 
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Paramagnetic solutions of the ionic Hubbard model at half-filling in dimensions D > 2 indi- 
cate that the band and the Mott insulator phases are separated by a metallic phase. We present 
zero-temperature dynamical mean-field theory solutions, which include antiferromagnetic long-range 
order, and show that the one-particle spectral functions always possess an energy gap and there- 
fore the system is insulating for all interaction strengths. The staggered charge density modulation 
coexists with antiferromagnetic long-range order of Neel type. 
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A bipartite lattice system of non-interacting electrons 
with one particle per site is a perfect gapless metal. Ap- 
plying an external alternating potential with a periodic- 
ity of twice the lattice constant a doubles the unit cell, 
thereby reduces the Brillouin zone (BZ), and opens a gap 
at the BZ boundary. Such a system is a band insulator 
with a charge density modulation with wavelength 2a; 
this system is also referred to as an ionic insulator^ On 
the other hand, switching on a local repulsive interaction 
between the electrons with opposite spins leads either to 
a paramagnetic Mott-Hubbard insulator with a correla- 
tion induced energy gap or to the spontaneous develop- 
ment of antiferromagnetic (AF) long range order. In the 
latter case it is the presence of AF order which doubles 
the unit cell and opens a gap at the BZ boundary for 
weak interactions and thereby creates a Slater insulator. 
In the strong interaction limit the electrons are local- 
ized with antiferromagnetically aligned spins forming a 
Mott-Heisenberg insulator. Experimental and theoretical 
investigations of metal-insulator transitions and transi- 
tions between different insulators continue as a challenge 
for condensed matter physics. 

The ionic Hubbard modeU^ incorporates both interac- 
tions and an external alternating potential and is there- 
fore well suited to study transitions between metallic 
or different insulating phases. This model was origi- 
nally used to study the neutral-ionic transition in or- 
ganic charge transfer salts 3 - or ferroelectric transitions in 
perovskite materials^ But the understanding of possible 
phase transitions in the ionic Hubbard model may prove 
important for other strongly correlated electron systems 
as well such as, for example, FeSi£ The physics of the 
ionic Hubbard model may even find a realization in op- 
tical lattices, if two laser beams of commensurate wave- 
lengths are superposed with properly tuned amplitudes. 6 

Extensive literature records exact, approximate, as 
well as numerical results for the ionic Hubbard model 
in one dimensioni 7 -^ It is agreed that at half-filling and 
in the interaction dominated regime the system is a para- 
magnetic Mott insulator whereas in the alternating po- 
tential dominated regime the system is an ionic band 



insulator. By now there is an emergent consensus that 
these two types of insulators are separated by yet another 
insulating phase with a non-zero bond-order parameter, 
which is the expectation value for a staggered component 
of the kinetic energy. 

More recently the ionic Hubbard model was also in- 
vestigated in higher dimensions within single-site or clus- 
ter dynamical mean-field theory (DMFT^iii^ii 3 . or by 
determinant quantum Monte-Carlo simulations?^ Using 
DMFT with iterated perturbation theory as the tool for 
solving the DMFT equations, Garg et al. determined the 
ground-state phase diagram of the ionic Hubbard model 
at half- filling for a semicircular density of states. 11 With 
the restriction to paramagnetic solutions an intermediate 
metallic phase was found separating the Mott and the 
band insulator. Within the same computational frame- 
work Craco et al. identified a coexisting phase between 
two insulators as well as discontinuous metal-insulator 
transitions*^ The discontinuous transitions were con- 
firmed in the two dimensional (2D) system by Kancharla 
et al.,— who used cluster DMFT combined with exact 
diagonalization and interpreted their data in favor of an 
intermediate bond ordered phase. 

While the DMFT work was restricted to paramagnetic 
solutions, finite temperature quantum Monte Carlo in 
2D simulations also probed AF correlations^ The pres- 
ence of the intermediate metallic phase between the band 
and the Mott insulator was confirmed. However, since in 
the 2D system at finite temperature long range antifer- 
romagnetism is prohibited^ the question how the possi- 
ble presence of the antiferromagnetic long-range order in 
higher dimensions or in 2D at zero temperature changes 
the phase diagram has remained open. 

Here we apply DMFT to the ionic Hubbard model at 
zero temperature allowing for spontaneous AF long-range 
order. We find that the ground state is always insulating, 
with a gap in the one-particle spectral function. There 
is a direct transition between the band insulator and the 
AF Mott insulator. Beyond a critical interaction strength 
the charge- and the spin-density modulations coexist. In 
this region the insulator has AF character. 



The ionic Hubbard model on a bipartite lattice is de- 
fined by the following Hamiltonian 
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where Aj = ±A/2 for i e A, B sublattices, respec- 
tively. The first term describes the kinetic energy of 
the electrons with spin a = ±1/2 for the hopping be- 
tween nearest-neighbor lattice sites i and j with ampli- 
tude t, the Hubbard interaction leads to an energy in- 
crease U > for the double occupancy of a site, and the 
last term contains a staggered potential with an energy 
difference A > between the A and B sublattices. The 
operators a ia and a) i(y obey standard fermionic anticom- 
mutation relations and = a) i(y a i(y . 

The Hamiltonian (p} is solved within DMFT at half- 
filling, i.e. for the chemical potential /j, = U/2, by a 
mapping to two non-equivalent impurity problems with 
the ionic energies ±A/2 . 16 i 17 ' 18 i 19 The impurity sites are 
coupled to two particle baths. Explicitly, we solve sepa- 
rately the two different single-impurity Anderson models 

H SIAM = ( e « _ M)»W + 22 V kot a {a ] aa C kacy + h.c) 

k 

(2) 
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where the hybridization matrix elements V kaa and the 
kinetic energies of the bath electrons t kaa are obtained 
self-consistently by additional DMFT equations^ Here 
we explicitly keep the spin and the site (a — A, B) depen- 
dences in order to allow selectively for spin and charge 
order. The one-particle impurity energy is e a — ±A/2 
for a = A or B, respectively. The single-impurity An- 
derson Hamiltonians are solved at zero temperature by 
the numerical renormalization group (NRG) method^ 
This method allows to obtain the spectral functions at 
and near the Fermi energy with high precision and can 
therefore accurately distinguish between metallic and in- 
sulating phases i^i 

The local (impurity) Green functions obtained from 
(|2"|) are expressed via the hybridization function r) acr (ui) 
and the self-energies Y> aa (w) as 



G aa (uj) 
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The hybridization functions describe the resonant broad- 
ening of the impurity energy levels due to the coupling 
to the particle baths and are given by 



Vkao 
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(4) 



The self-energies capture the interaction induced corre- 
lation effects on the impurity sites. 

Within DMFT the hybridization functions are subject 
to self-consistency conditions, which involve the density 




0.8 
0.4 



(E) 



U=2.5 



-3-2-10123 
Energy, CO 

FIG. 1: (color online) Spectral functions for the ionic Hub- 
bard model at A = 0.5 and different interactions U in the 
paramagnetic limit. Solid and dashed lines correspond to the 
A and B sublattices. From top to bottom: (A) band insulator; 
(B), (C) correlated metal; (D), (F) correlated Mott insulators. 
Note the different scales on the axis in particular horizontal 
ones. 



of states (DOS) for a given lattice structure ^ 7 ! 18 ' 19 In the 
following we adopt the semicircular DOS corresponding 
to the Bethe lattice ! 19 ! 22 ' 23 The hybridization functions 
are then simply related to the local Green functions © 
through 
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where a = B , A if a — A, B and a = —a, respectively. 
IF = 1 is the bandwidth, which sets the energy unit. 
We emphasize that the use of the semicircular DOS has 
merely technical reasons because this choice simplifies the 
DMFT equations. The obtained results remain qualita- 
tively similar for any particle-hole symmetric DOS, which 
represents a bipartite lattice in dimensions D > 2 i 22 i 23 

From the self-consistent DMFT solution of the ionic 
Hubbard model we determine (i) the local one-particle 
spectral function A aa {uj) = — ImG Q , IT (a;)/7r, (ii) the 
AF order parameter (staggered magnetization) toaf = 
(n^4| — riAi) = — (tib] — "-si), and (hi) the charge den- 
sity wave amplitude (ionicity) tticdw = { n A-\ + rc-Aj. — 
riB] — tibi), where (• • •) denotes ground-state expecta- 
tion values. A metal is distinguished from an insula- 
tor by a finite spectral function at the Fermi level, i.e. 

^(o) = E w A«r(o)>a 

Without allowing for long-range AF order in the self- 
consistent DMFT solution the results of iterated pertur- 
bation theory suggested that the band and the Mott in- 
sulator are separated by a metallic phase. 11 ! 12 Here, we 
confirm this conclusion by using NRG for solving the 
DMFT equations. In Fig.[T]we show spectral functions in 
the paramagnetic regime at fixed A = 0.5 starting from 
a band insulator at small interactions (cf. top panel with 
U = 0.5). By increasing U we find a continuous tran- 
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FIG. 2: (color online) Spectral functions at the Fermi energy 
versus U for the ionic Hubbard model at A = 0.5 in the para- 
magnetic limit. Solid data points (black and red) are obtained 
from DMFT iterations which start from an initial metallic in- 
put whereas open points (blue and green) are obtained from 
an insulating input. The band insulator to metal transition 
is continuous, but the metal to Mott insulator transition is 
hysteretic. Inset: charge density wave amplitude mcDW vs. 
U. 
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FIG. 3: (color online) Spectral functions of the ionic Hubbard 
model for A = 0.5 and different interactions U allowing for 
AF order. Solid and dashed lines correspond to different sub- 
lattices, the color code distinguishes the spins a = ±1/2, as 
indicated in the inset. (A) Paramagnetic band insulator with 
meow > and jtiaf = 0. (B) - (E) Antiferromagnetic ionic 
insulator with mcDW > and ttzaf > 0. Note different scales 
on vertical axis. In all cases the spectra are gapped. 



sition to a metallic solution with finite spectral weight 
at the Fermi energy (for U — 1 and 1.5 in Fig. QJ3 and 
C). This metallic phase coexists with long-range charge 
order, i.e. tocdw > 0. Larger interaction strengths ho- 
mogenize the system and the charge density wave ampli- 
tude continuously decreases (cf. the results for U = 1.5 
in Fig. [Tp and the inset in Fig. [2]). A further increase of 
U leads to a Mott-Hubbard type metal-insulator transi- 
tion with hysteretic behavio n 21 i 24 In the examples shown 
in Fig. Q] for U — 2 and 2.5 we find Mott insulators with 
correlation induced spectral gaps. In the Mott insulator 
the charge density wave is very small but remains finite. 

The two transitions are separated in the phase dia- 
gram and imply the existence of two critical interaction 
strengthsj 11 ' 12 The transition from the metal to the Mott 
insulator resembles the one, which is found within DMFT 
applied to the paramagnetic Hubbard model in the ab- 
sence of a staggered potentiah 21 i 24 At zero temperature 
this transition is continuous, though hysteretic behav- 
ior in the iterative solution is encountered; it occurs at 
U c w 1.45 for a semicircular DOS. 

We note that although tocdw becomes vanishingly 
small for large U (see the inset in Fig. [2]) it is expected 
to remain finite for all interaction strengths as long as 
A > 0. This expectation relies on the general argument 
that the symmetry of the ground state cannot be higher 
than the symmetry of the Hamiltonian itself. 

It is understood, however, that the paramagnetic solu- 
tion of the Hubbard model at half-filling is not the generic 
case. Antiferromagnetism is likely to occur unless it is 
prohibited e.g. by strong frustration effects. In the next 
step we therefore discuss the ionic Hubbard model at 



half-filling on a bipartite lattice allowing for long-range 
antiferromagnetism. 

Both, staggered charge order, induced by the alter- 
nating potential A, and spontaneous staggered AF or- 
der may give rise to a gapped spectrum and insulating 
behavior of the lattice system at half-filling. If both 
staggered orders develop simultaneously, the above men- 
tioned metallic phase in between the ionic band and the 
Mott insulator may be insulating due to an AF induced 
energy gap. Indeed, as the examples in Fig. [3] illustrate, 
the spectral functions are gapped around the Fermi en- 
ergy in the entire U range. At weak U the system pos- 
sesses charge order only; staggered magnetization is zero 
(cf. Fig. By increasing U tucdw is reduced and 

AF correlations develop. The local interaction U reduces 
the amount of double occupancy - needed to maintain a 
finite mcDW at half-filling. On the other hand, virtual 
hopping processes induce an effective exchange interac- 
tion, which favors antiferromagnetism. Beyond a critical 
(7(A) the system acquires Neel order (cf. Fig. [3]B - E) 
accompanied by a strong reduction of the CDW ampli- 
tude. 

In Fig. 0] we present how the order parameters tocdw 
and toaf vary with the interaction U for different ionic 
potentials A. As discussed above, due to the suppression 
of double occupancies by the on-site repulsion, the charge 
density wave order parameter is reduced but tocdw > 
for all U as inferred from symmetry arguments. Whereas 
at A = Neel order appears at infinitesimally small U, 
at finite A the interaction has to exceed a finite critical 
value U C (A) for the onset of antiferromagnetism. This 
quantum phase transition at U = U C (A) is continuous 
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FIG. 4: (color online) U dependence of spin- and charge- 
density wave order parameters for different ionic potentials 
A. 

in contrast to the paramagnetic-antiferromagnetic tran- 
sition for models with frustration.— At very large U the 
AF order parameter saturates at its maximum value. 
In summary, our numerical solution of the DMFT 



equations for the ionic Hubbard model provides evidence 
for the existence of a critical interaction strength for the 
transition from a weakly correlated band insulator to a 
Mott insulator with coexisting charge and staggered spin 
order. A gap in the one-particle spectrum persists in all 
parameter regimes and thus implies the absence of an 
intervening metallic phase. We emphasize that our re- 
sults do not contradict the findings in Refs. [Tummi 

because in those works the AF long-range order was ei- 
ther excluded by the choice of the method or by the 
low dimensionality of the system and finite temperatures. 
Transitions between a correlated metal and a band insu- 
lator are also found in other models, such as the Hub- 
bard model with binary alloy disorde r 26 i 27 or the bilayer 
Hubbard model with interlayer hopping.— As we have 
demonstrated here, the occurrence of spontaneous stag- 
gered long-range order can significantly change the na- 
ture or even the existence of a transition to a metallic 
phase. 
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of the Deutsche Forschungsgemeinschaft (DFG). The 
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